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ABSTRACT 

Oil consumption was measured in a modem production diesel engine using tritium as a 
radiotracer. The measurements were made primarily at two speeds and one load using first 
a single-grade lubricant and then a multi-grade lubricant. 

These values were then compared to oil flow rates up/down the liner which were based 
on film thickness traces of a sister engine under the same loads and speeds. The traces 
were obtained using the laser-fluorescence technique. For the most part, it was discovered 
that there does not seem to exist a correlation between these flow rates and oil 
consumption. However, the traces do reveal that the crown land is dry on all four strokes 
and thus does not contribute to the engine's oil consumption. 

A larger data base is necessary in order to accurately compare oil consumption to the 
film traces. This is currently in progress as of this writing. 

Thesis Supervisor: David P. Hoult 

Title: Senior Research Associate 
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CHAPTER 1 BACKGROUND AND INTRODUCTION 

Due to the EPA particulate emission standards for 1991 and 1994, it has become 
necessary to learn more about oil consumption and the driving forces behind it. 
Subsequently, an accurate determination of these driving forces could lead to designs of 
pistons, rings, and lubricants which improve both oil consumption and particulate 
emission. 

This paper, then, has two purposes. First, initiate an industry proven method for 
measuring oil consumption, and second, compare these values to the average oil flow rates 
up the liner of a sister engine operating at the same loads and speeds. The oil consumption 
measurement technique, along with engine set-up and experimental procedures, will be 
explained further. The oil flow rates, however, have been calculated from oil film 
thickness data taken from, and explained in, separate papers [1,2 respectively]. 
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CHAPTER 2 DESCRIPTION OF MEASURING TECHNIQUE 

The technique used is a radiotracer method which uses tritium as its radiotracer. This 
method has been widely used for assessing oil consumption for over two decades. The 
method is fast, accurate, and economical relative to non-radioactive methods. This overall 
practicality and effectiveness is a major result of the advances of liquid scintillation 
counting. 

2.1 BASIC METHOD 

The following is a brief description of the basic method while a more detailed 
description of the general technique is found in Reference [3]. However, for ease of 
reference, the derivation of the oil consump*ion calculation equation, upon which this 
method is based, is included as Appendix A. 

The basic procedure consists of running an engine with radioactive lubricating oil with 
a known radioactive level of disintegrations per minute (dpm). The total amount of 
hydrogen passing through the engine, via fuel flow and air flow, is measured and 
calculated. Through combustion, this hydrogen, along with the hydrogen/tritium from the 
oil, leads to the formation of water. A sample of this water is collected by condensing a 
sample of the exhaust gas. (Only a sample is needed since total measurements of the flow 
rates into the engine are obtained) Subsequently, the radiation concentrations of both the 
condensed water and the oil are determined using a liquid scintillation counter. These 
values are then applied to the oil consumption equation to arrive at the engine's total oil 
consumption. 

2.2 SOOT SAMPLE COLLECTION 

As seen in Appendix A. the general technique is based on complete combustion of the 
exhaust sample. For this paper, however, the exhaust sample is filtered and then 
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CHAPTER 3 EQUIPMENT AND SET-UP 

3.1 ENGINE 

The engine used in this experiment is the KUBOTA EA300N, a single cylinder, 4- 
stroke, IDI diesel. This engine is used primarily for remote power generation. A detailed 
description of the engine's geometry and performance is included in Appendix B. 

This engine was chosen only because the film thickness measurements used in this 
paper were obtained from an engine of the same make and model [1]. 

3.2 RADIOACTIVE LUBRICATING OILS 

Radioactive oil for these experiments is purchased commercially. The oil received has 
been through a tritiation process which results in approximately 250 mCu in one ml of oil. 
Subsequently, this oil is divided up and thoroughly mixed with Pennzoil SAE-30 (single- 
grade) and Pennzoil 15W-40 (multi-grade). The tritiation process is explained further in 
Appendix B . 

These oils are used primarily for the same reason as above for the engine. In addition, 
they should provide a good comparison of single- and multi-grade oils since the only 
difference is the addition of a polymer VI improver to the 15W-40.[3] 

3.3 ENGINE INSTRUMENTATION 

The engine is fitted with thermocouples ji order to monitor operating temperatures (Fig. 
1). The thermocouples measure the inlet air temperature, the coolant temperature, the oil 
sump temperature, the exhaust temperature, and the liner temperatures. The latter 
thermocouples coincide with the top ring at TDC, mid-stroke, and BDC. A laminar air 
flow element is also used to measure the air flow while a burette system is used to measure 
the volume flow rate of the fuel. A more detailed listing of the measuring equipment is 
included in Appendix B. 



12 



condensed instead of being completely combusted and then condensed. This results in an 
activity concentration (dpm/ml) of the condensed water and a total activity (dpm) of the soot 
sample. 

In order for this method to be valid, two assumptions are made. First, the soot 
collected, if burned, would produce a negligible amount of water. Secondly, all of the 
water vapor in the sample is being condensed. This latter assumption is necessary since the 
consumption calculation requires radioactive concentrations vice total activities. Therefore, 
the soot sample total activity is divided by the amount of water collected and then added to 
the activity concentration of the water sample. A more in depth description of the 
procedures is included in Appendix E. 

2.3 LIQUID SCINTILLATION COUNTING [41 

All radioactive levels are determined by a method known as liquid scintillation counting. 
This technique is highly sensitive and greatly responsible for the increased use of weak beta 
emitters such as carbon-14 and tritium. The basic technique, as with other methods of 
measuring radioactivity, relies on the interaction of nuclear radiation with matter. In this 
case, a luminescent material, or scintillator, is excited by radioactive disintegrations. 

A scintillator, also referred to as either a phosphor or sometimes a fluor, is defined as a 
material which emits a brief pulse of fluorescent light when interacting with a high-energy 
particle or quantum. Therefore, the radioactive samples used in this experiment are mixed 
with a scintillator. The subsequent flashes of light, or scintillations, are detected by the 
counter by use of a photomultiplier. This in turn, over a set period of time, allows for the 
determination of the average disintegrations per minute (dpm) of each sample. 
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3.4 DESCRIPTION OF PISTON AND RINGS [21 

The piston and rings (Fig. 2) are representative of modem engine design. The piston is 
made of aluminum and accommodate three piston rings. It has round lands and a barrel 
shaped skirt. Behind and below the oil control ring, there are oil flow relief holes that 
allow the circulation of oil up the piston skirt on the liner side and back down to the sump 
inside the piston. 

The three piston rings consist of a top ring, a scraper ring, and an oil control ring. The 
top ring is a half-keystone design with a 7° keystone angle. It is made out of cast iron with 
a chromium strip embedded into its face to improve wear resistance. The resulting face 
profile is asymmetric. The scraper ring has a tapered face and an undercut on its face side 
to provide the correct tilt angle when compiessed to its working diameter. It is also made 
of cast iron. The oil control ring is a one price type made of cast iron also, and, like the top 
ring, its rails have a chromium strip embedded into them. The oil control ring also has 
slots between the rails to allow the flow of oil through the ring to the holes in the piston 
mentioned above. 

3.5 SAMPLING SYSTEM 

Figure 3 is a schematic of the sampling system in use. As shown, the sampling system 
is connected to an exhaust tank. This tank is used to help eliminate the exhaust pulse that 
exists when operating a single cylinder engine. This entire system, consisting of 
inexpensive components, is placed on a small lab cart which allows for portability and for 
an easy connection to a different engine. 
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3.6 OIL FILM THICKNESS MEASUREMENTS 

The film thickness traces referred to are obtained from reference [1] which uses a laser- 
fluorescence technique. This basic technique consists of mixing the lubricant with a 
fluorescent dye and then, during operation, focusing a HeCd (blue) laser through a quartz 
window installed in the liner. The resulting fluoresced (green) light is then collected and 
converted to a voltage signal using a photomultiplier [2], Lastly, the voltage signal is 
converted into microns using a calibration coeffecient (unique for each trace). This 
coeffecient is based on the voltage readings obtained when etch marks of known depths on 
the piston skirt pass the quartz window.[l] 

More in depth descriptions of this technique can be found in references [1] and [2]. 
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CHAPTER 4 DESCRIPTION OF EXPERIMENTS 

4.1 PROCEDURES 

At each given load and speed, the engine is allowed to stabilize. Subsequently, flow 
rates are measured and samples are taken in accordance with the step-by-step procedures in 
Appendix C. The samples are then prepared in accordance with Appendix D in order for 
the scintillation counter to be effective. 

Following sample preparation, the samples are placed in the scintillation counter to 
determine their individual radioactive levels. These values, along with the measured flow 
rates, are then used to calculate oil consumption. Actual formulas and an example of the 
spreadsheets in use are included in Appendix E. 

After collecting data using the first oil, the engine is flushed to minimize contamination 
of the next test oil. A flow chart of these basic procedures is illustrated in Figure 4. 

4.2 ENGINE OPERATING CONDITIONS 

The primary loading conditions for each oil consisted of full load at 1500 rpm and 
full load at 3000 rpm (approximately 9.5 ft-lbs for each). As before, these conditions were 
used in order to duplicate the conditions under which the Film thickness data was taken. 
Secondary loading conditions consisted of 2000 rpm and 2500 rpm at the same torque. 
These latter conditions were used mainly for graphic continuity between low and high 
speeds. 
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CHAPTER 5 GRAPHIC RESULTS 

5.1 OIL CONSUMPTION MEASUREMENTS 

Oil consumption measurements for the desired load and speeds using Pennzoil's SAE- 
30 and 15W-40 are illustrated in figures 5 and 6, respectively. Figure 7 is then an overlay 
of the previous two graphs with figure 8 being another comparison but in the form of a bar 
chart. 

The volatility of a new oil (SAE-30) is also examined. The results, included as figure 9, 
illustrate how oil consumption rates decrease and then stabilize with time (for a constant 
load and speed). Volatility information is desirable since the film thicknesses were 
obtained from the sister Kubota using new lubricating oils. 

5.2 OIL FILM THICKNESSES 

Figures 10, 11, 12, and 13 illustrate a subset of the average film thickness traces for ten 
compression, expansion, exhaust, and intake strokes, respectively [1]. For these traces, 
the engine was using SAE-30 and was running at full load at 1500 rpm. The traces for 
3000 rpm and for the different oil are similar to these and are therefore not included. 
Interpretation of this data, which can be confusing to the first-time observer, is described in 
length in reference [2], 
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CHAPTER 6 ANALYSIS ANI> DISCUSSION 
6.1 HYPOTHESES 

There exists three prevalent hypotheses on lubricating oil consumption, two of which 
deal with oil being drawn in and around the top ring at or near the end of the expansion 
stroke. At this point, the high pressure below the top ring forces the oil around the top ring 
where it either (1 ) accumulates on the crown land from where it is later burned off or (2) is 
"blown" up into the combustion chamber and consumed. The third hypothesis is based on 
oil being transported by the top ring face. Subsequently, some of this oil is deposited onto 
the liner at or near top ring reversal (TRR) where it is later consumed. 

For each hypothesis, there exists a corresponding volume. These volumes are 
calculated by first finding their respective areas which are illustrated in figure 14. The areas 
are then swept around the circumference of the bore to determine the volumes. Table 1 lists 
the areas and the volumes calculated. 



SPACE 


AREA (mm 2 ) 


VOLUME (min^) 


A) Crown land crevice 


1.33 


313.0 


B) Ring/groove crevice 
(top and bottom) 


0.3 


70.7 


C) Ring/groove crevice 
(backside of ring) 


1.02 


241.0 


D) Between ring and liner 


.00386 


0.91 



Table 1. Areas A^olumes of possible oil consumption mechanisms 
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These volumes can in turn be related to an oil flow rate, in g/hr, based on engine 
operating speed and oil density (Refer back to Appendix E). Unfortunately, volumes A, B, 
and C above result in rates on the order of thousands of grams per hour with volume D 
being on the order of 50 grams per hour. These values are extraordinarily high and do 
not scale with oil consumption (figures 5 and 6). Thus they cannot offer much guidance in 
determining oil consumption mechanisms. Actual volumes necessary to resemble the oil 
consumption values are on the order of 0.01% of volumes A and C, 0.05% of volume B, 
and just 3.0% of volume D. 

6.2 ANALYSIS AND DISCUSSION 

As discussed in Chapter 5 of reference [2], the film thickness data does not usually 
allow for a determination of where the oil lies (i.e. on the piston, on the liner, or on both). 
However, as figures 10 through 13 illustrate, the oil film thickness before the piston arrives 
(-4mm to 0mm for the upstrokes) is relatively the same as when the crown land is passing 
(0mm to 9.5mm for the upstrokes). Therefore, it is apparent that the piston crown land is 
virtually oil free and thus, does not contribute to oil consumption. It is important to 
mention that this film trace characteristic is also shared by the other traces for the higher 
speed and with the different oil. 

Based on the above line of reasoning, the first hypothesis can be easily refuted. The 
second one, however, can neither be directly refuted nor explicitly proved by this method 
of film measurement. This is because the hypothesis addresses oil which enters the 
combustion chamber by way of neither the piston nor the liner. However, further data 
analysis of the film traces does allow for some broad conclusions to be made concerning 
both this and the final hypothesis. 

This subsequent data analysis first requires the use of existing computer programs 
which have been modified for the Kubota engine [5]. These programs calculate the volume 
differences (mm-Vcycle) of various crank angle regions for both the power and the gas 
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exchange strokes. Only five regions, which are listed below in Table 2, are analyzed in 
this paper. The methods and reasons for computing these particular regions, along with 
copies of the programs, are included in Appendix F. 



REGION 


PISTON AREA 


I 


Crown land 


II 


Top Ring 


m 


1 ring width below top ring 


IV 


Second Land 


V 


Second (Scraper) Ring 



Table 2. Piston Regions of Film Thickness 



The title of each region in the table explains which part of the piston is passing the 
quartz window at that instant. In addition, the regions are labeled on an expanded version 
of the aforementioned compression stroke. (Figures 15 and 16) 

Secondly, in order to compare a regions volume difference to the measured oil 
consumption, the former is converted to an oil flow rate up (or down) the liner in g/hr. 
This is done based on the volume difference, the speed, and the oil's density. (Refer back 
to Appendix E for calculations) 

Tables 3 and 4 and Tables 5 and 6 tabulate the data collected using SAE-30 and 15W- 
40, respectively. A positive value indicates that the upstroke volume is greater than the 
downstroke and vice versa. Subsequently, figures 17 through 22 are included to illustrate 
these various flow rates of the two oils for the different speeds. On each figure, the 
respective oil consumption is included in order to make relative comparisons. 
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SAE-30 at 1500rpm 


PWR EXCHANGE 


GAS EXCHANGE 


TOTAL (pwr+gas) 




Vol Diff 


Oil Flow 


Vol Diff 


Oil Flow 


Vol Diff 


Oil Flow 


Region 


I 


0.04 


1.51 


0.00 


0.00 


0.04 


1.51 


Region 


II 


0.01 


0.38 


-0.08 


-3.03 


-0.07 


-2.65 


Region 


m 


1.11 


42.01 


-0.23 


-8.70 


0.88 


33.30 


Region 


IV 


2.30 


87.04 


-0.52 


-19.68 


1.78 


67.36 


Region 


V 


-0.02 


-0.64 


-0.11 


-4.16 


-0.13 


-4.81 



Table 3. Oil Flow Rates up the Liner (SAE-30 at 1500 rpm) 



SAE-30 at 3000rpm 


PWR EXCHANGE 


GAS EXCHANGE 


TOTAL (pwr+gas) 




Vol Diff 


Oil Flow 


Vol Diff 


Oil Flow 


Vol Diff 


Oil Flow 


Region 


I 


-0.34 


-12.87 


-0.12 


-4.54 


-0.46 


-17.41 


Region 


II 


0.08 


3.03 


-0.21 


-7.95 


-0.13 


-4.92 


Region 


III 


0.74 


28.01 


-0.23 


i 

OO 


0.51 


19.3 


Region 


IV 


1.43 


54.12 


-0.45 


-17.03 


0.98 


37.09 


Region 


V 


0.32 


12.11 


-0.19 


-7.19 


0.13 


4.92 



Table 4 Oil Flow Rates up the Liner (SAE-30 at 3000 rpm) 



22 



15W-40 at 1 500rpm 


PWR EXCHANGE 


GAS EXCHANGE 


TOTAL (pwr+gas) 




Vol Diff 


Oil Flow 


Vol Diff 


Oil Flow 


Vol Diff 


Oil Flow 


Region 


I 


-0.04 


-1.51 


0.08 


3.03 


0.04 


1.51 


Region 


II 


0.03 


1.14 


-0.04 


-1.51 


-0.01 


-0.38 


Region 


in 


0.63 


23.84 


0.05 


1.89 


0.68 


25.73 


Region 


IV 


1.17 


44.28 


-0.08 


-3.03 


1.09 


41.25 


Region 


V 


0.16 


6.06 


-0.01 


-0.38 


0.15 


5.68 



Table 5. Oil Flow Rates up the Liner (15W-40 at 1500 rpm) 



1 5W-40 at 3000ipm 


PWR EXCHANGE 


GAS EXCHANGE 


TOTAL (pwr+gas) 




Vol Diff 


Oil Flow 


Vol Diff 


Oil Flow 


Vol Diff 


Oil Flow 


Region 


I 


-0.22 


-8.33 


-0.16 


-6.06 


-0.38 


-14.38 


Region 


II 


0.10 


3.78 


-0.20 


-7.57 


-0.10 


-3.78 


Region 


in 


0.46 


17.41 


-0.38 


-14.38 


0.08 


3.03 


Region 


IV 


0.91 


34.44 


-0.66 


-24.98 


0.25 


9.46 


Region 


V 


0.09 


3.41 


-0.35 


-13.25 


-0.26 


-9.84 



Table 6. Oil Flow Rates up the Liner (15W-40 at 3000 rpm) 



Upon examination of these Figures, it is seen that a majority of the regions show 
absolutely no correlation to oil consumption. However, a slight correlation does exist for 
Region II, the top ring, but only during the power exchange strokes. Therefore, a closer 
look is taken at non-dimensional ized values for this case and for the oil consumption. This 
is accomplished by dividing each rate by its corresponding 1 500 
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rpm value, (i.e. all of the 1500 rpm values will equal one) The results are illustrated in 
figure 23 for SAE-30 and figure 24 for 15W-40, however, it is seen that the oil flow rates 
do not duplicate the trends of the oil consumption in either case. 



6.3 PROBABLE ERRORS IN ANALYSIS TECHNIQUE 

When comparing oil flow rates to oil consumptions, if in fact film thickness in this case 
can even be related to oil consumption, only trends can be looked for. An exact matching 
of values would be virtually impossible for the following reasons: 

a) Oil consumption measurements vary within an engine, therefore they will 
vary even more between separate, but identical, engines. 

b) The oil film data is obtained white using new, and thus more volatile, 
lubricating oil. (Refer back to figure 9) 

c) As previously mentioned, the calibration coefficients for each film trace are 
determined by way of etch marks on the skirt [1], However, there may 
exist a temperature gradient between the skin and the rings and lands. 

This gradient would result in different calibrations for these latter 
regions since the fluorescence of the oil is inversely proportional to its 
temperature [6]. 

d) The oil flow relief holes, described in sect. 3.4, is probably what attributes 
to the absence of correlation between oil consumption and film thickness. 
Figure 25 illustrates that, over the length of the piston, the power strokes 
tend to have more oil on the way up then on the way down and vice versa 
for the gas exchange strokes. TTiese results imply that: 

1 ) after the oil control ring passes the window on the compression 
stroke and before it passes the window on the expansion stroke, oil 
is returned to the sump via the flow holes. 

2) after the oil control ring passes the window on the exhaust 
stroke and before it passes the window on the intake stroke, oil 
is supplied to the sump via the flow holes. 
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6.4 MISCELLANEOUS ANALYSIS 

The temperature of the upper portion of the liner during operation is plotted in figure 26 
for the RPMs shown. These temperatures are then used to determine the high and low 
shear viscosities (figure 27) and the surface tension (figure 28) of the oils in 
use [7]. Lastly, the oil consumptions versus viscosities are plotted in figure 29. 
Unfortunately, no trends or relationships between the two are observed. 
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CHAPTER 7 CONCLUSIONS AND RECOMMENDATIONS 

The primary conclusions and recommendations are as follows: 

1) The crown land is virtually dry and therefore does not contribute to oil 
consumption. 

2) The top ring is the only region whose oil film thickness characteristics remotely 
resemble oil consumption. However, before any conclusions can be formulated, more 
comparisons need to be made between this region's flow rate and oil consumption. 

3) The concept of oil being "blown ' around the top ring and into the combustion 
chamber is still a viable solution, however, the top ring region might be an additional 
driving force behind oil consumption. 

4) The oil flow relief holes provide a significant unknown which might make this 
type of flow rate analysis for this particular set-up invalid. A recommended set-up would 
consist of obtaining film thickness measurements through a window placed such that the oil 
control ring never passes over it. 
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APPENDIX A DERIVATION OF OIL CONSUMPTION 
CALCULATION 131 



The following derivation is based on the composition of the combustion water. This 
water comes from the water in the inlet air and from its formation via the combustion 
process of the fuel and of the oil. If radioactive hydrogen is only in the oil, then: 



disintegrations per minute (dpm) = Lube oil consumption dpm (a.l) 
Wt of sample of combustion H2O H2O from fuel comb. + H2O 

from oil comb. + H2O from air 



The following symbols are used: 

R 0 = Wt rate of oil consumption, g/hr 
R f = Wt rate of fuel consumption, g/hr 
R a = Wt rate of air consumption, g/hr 

K 0 = Water fonnation constant of oil, g water formed/g oil burned 
Kf = Water fonnation constant of fuel, g water formed/g fuel burned 
h s = Specific humidity, g H 20 /g air 

SA ov = Specific Activity for H2O sample, dpm/g of sample 
SA 0 = Specific Activity for lube oil, dpm/g 
t = time, hr 

These values, when substituted into equation, a.l. result in: 

= Rq( (S Aq )/[ Rq( Ky + Rf( Kf + R a thg] 
Subsequent grouping of the R 0 terms results in: 

R 0 = [R f K f + R a h s ]/[(SA 0 /SA CW )-K 0 ] 
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The water formation constants, K 0 and Kf, are detennined based on their hydrogen to 
carbon ratios (H/C). (See Appendix E) 

The H/C ratios for the Pennzoil oils are obtained via Dr. Tayeb Benchaita of Pennzoil 
Products. The H/C for the fuel is based on a standard light diesel fuel. Actual values and 
sample calculations are included in Appendix E. 
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APPENDIX B DETAILED EQUIPMENT DESCRIPTION 
B.l KUBOTA EA300N CHARACTERISTICS 

The Kubota EA300N's detailed geometry is presented in the following table. 



Make 


Kubota 


Model 


EA300N 


Type 


Horizontal, 4-stroke, IDI Diesel 


Number of cylinders 


1 (one) 


Bore and Stroke 


75 x 70 mm (2.95 x 2.76 in.) 


Displacement 


0.309 liters (18.86 cu. in. ) 


Compression Ratio 


23:1 


Cooling System 


Water cooled - natural convection 


Lubrication System 


Trochoidal pump (no oil filter) 


Rated Brake Horsepower 


4.48 kW (6hp) @ 3000 RPM 


Maximum Torque 


15.2 Nm (11.2 ft-lbf) 


Typical Application 


Remote Power Generation 



Table B.l Engine Characteristics 



A unique characteristic of this engine is that it does not have a paper filter. Instead, it 
has a strainer screen and a magnet to keep large pieces of debris and metal shavings from 
entering the oil pump. Unfortunately, this allows the dirt and grit that is introduced into the 
oil during operation to remain in the oil. Therefore, the engine oil becomes black very 
quickly. This "blackness" requires the oil samples to be diluted and centrifuged (see 
Appendix D) in order for the scintillation counting technique to be of use. 
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B.2 RADIOACTIVE OIL [61 

The tritiated oil received has been through a reduction process using tritium gas. This 
gas, plus one ml of the neutral base oil, was mixed with forty mg of a platinum black 
catalyst and stirred overnight. Afterwards, any labiles were removed using 50% benzene 
methanol. The final result is the addition of tritium atoms to the hydrocarbon molecules in 
the oil. This process is referred to as either a hydrogenation or tritiation process. 

B.3 ENGINE INSTRUMENTATION 

The following is a listing of the instrumentation used while collecting data for these 
experiments. (Refer back to Fig. 1) 



COMPONENT: 


Make, Model #: 


Liner Thermocouples (3) 


Omega, J-type 


Liner Thermocouple Readout 


Omega, Model #2168A 


Other Thermocouples 


Omega, K-type 


Thermocouple Readout 


Omega, Model #199 


Hall-effect Speed Sensor 


Minarik, Visi-Tach 


Load Cell 


Eaton, Model #3169 


Laminar Air Flow Element 


Meriam Instr., Model 50MW20-2 


Liquid Scintillation Counter 


Packard, 2000 CA TRI-CARB 



Table B.2 Measuring Equipment and Instrumentation 
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APPENDIX C RADIOTRACER PROCEDURES 
C.l MIXING AND "BREAKING-IN" THE OIL 

Non-radioactive oil should be mixed with radioactive oil until the concentration level is 
approximately six microcuries per gram (6 jiCu/g). If a completely new batch of oil is 
used, the engine should be operated for around 5 hours before any data samples are 
collected. This break-in period allows any light ends to be burned off. (Refer back to 
figure 9) 

C.2 SAMPLING PROCESS 

Before samples are taken, operate the engine continuously at the given load and speed to 
allow temperatures to stabilize. Between speed changes, allow approximately 15 minute 
for stabilization. After stabilization, run the vacuum pump for a short period of time in 
order to draw exhaust gas through the sampling system. This is done (1) to ensure that the 
engine and sampling system are making good water samples and (2) to ensure that any 
water residue from the last data run is emptied from the condenser coils. The following is a 
step-by-step description of the actual sampling process. 

1 ) Replace water collection bottle with a clean one 

2) Replace exhaust filter with a new one 

3) Turn on the vacuum pump to initiate the sampling process 

4) Measure and record the volume flow rate of fuel using the installed burette 
and a tuner 

5) Record the pressure difference across the laminar flow element using the 
installed inclined manometer 

6) Record operating temperatures 

7) After approximately ten minutes, secure from sampling by turning off the 
vacuum pump. This time limit allows for. 
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a) collection of a sufficient amount of condensed water and 

b) a sufficient number of piston ring revolutions and thus a better 
average oil consumption 

8) Remove both the filter and the water collecting bottle and replace with clean 
ones (see Appendix D for preparation) 

9) Stabilize at the next load and speed and repeat from step 3 on. 

10) After running at the desired load and speeds, secure the engine and 

drain approximately 50 ml of lubricating oil followed by draining one ml of oil 
into a separate test tube (see Appendix D for preparation) 

1 1 ) Return the 50 ml of drained oil to the sump 

12) Record the barometric pressure and relative humidity 

All data is recorded on available data sheets in the units indicated on the sheet. Using 
these units allows the values to be inputted directly into the present spreadsheets (Appendix 
E). Copies of the data sheets in use are included as the next two pages. 
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EXPERIMENTAL D AT A SHEET 

U NIFORM INPUT : 

DATE 



Oil brand 



Oil Activity 

(dpm/.5 ml kero/oil so In) 



Barometric pressure 

(mm mercury ) 



Air temperature 
(cel si us ) 



Humidity coef f . 



Relative Humidity 
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Load 

(lbs) 


Load 

(lbs) 


Speed 

( rpm ) 


Speed 

( rpm) 


Fuel rate 

( ml /sec ) 


Fuel rate 

(ml/ sec ) 


Air press, diff. 
across element - 
(inches wtr) 


Air press . diff. 
across element - 
( inches wtr ) 


amount of water 

collected 

(ml ) 


amount of water 

collected 

(ml ) 


activ of soot 

sample 

(dpm ) 


activ of soot 

sample 

( dpm ) 


activ cone, of 

wtr sample 

( dp m/ml ) or ( dpm/g) 


activ cone . of 

wtr sample 

( d pm/ ml ) or ( dpm/g ) 


TEMPERATURES: 


TEMPERATURES: 


exhaust 

air inlet 

oil 

coolant 


exhaust 

air inlet 

oil 

coolant 


liner: 

upper 

middle 

lower 


liner: 

upper 

middle 

lower 
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APPENDIX D SAMPLE PREPARATION [71 
D.l WATER SAMPLES 

1) Filter the condensed water into a clean test tube 

2) Add and mix a sufficient amount of ammonium hydroxide solution to the 
sample to bring the PH up to 7 (neutral). [7,8] 

3) Centrifuge the sample for approximately ten minutes 

4) Filter again if necessary 

5) Pipette 10 ml of fluor solution in:o a counting bottle 

6) Pipette 1 ml of the water sample into the counting bottle (note - more or 
less water may be counted but the spreadsheet is set up for 1 ml —see 
Appendix E) 

7) Replace top, mix, and number the sample 

D.2 SOOT SAMPLES 

1) Place the filter in a tube containing approximately 30 ml of fluor solution and 
shake vigorously to mix well 

2) Allow two hours for the fluor solution to dissolve the soluble hydrocarbons 
and then shake again 

3) Filter the solution into a clean test tube. Repeat if necessary 

4) Pipette 10 ml of the fluor solution into three counting bottles each 

5) Divide the filtered solution among the bottles 

6) Replace the tops, mix, and number the samples 

D.3 OIL SAMPLES 

1) Pipette 10 ml of kerosene into a clean test tube 

2) Pipette .2 ml of the oil sample into the kerosene. Subsequently, flush 
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(he pipette with the oil/kerosene solution until all visible oil 

residue is 

removed from the pipette. Place pipette into the solution 

3) Mix well and centrifuge for approximately ten minutes 

4) Pipette 10 ml of the fluor solution into a counting bottle 

5) Pipette .5 ml of the oil/kerosene solution into the counting bottle (note - 
as before, more or less may be counted but the spreadsheet is set up for .5 
ml — see Appendix E) 

6) Replace top, mix, and number the sample 

All samples are taken to the Radiation Protection Office so that their activity levels may 
be counted using a liquid scintillation counter. 
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APPENDIX e formulas and example calculations 

E.l OIL CONSUMPTION 

E.la MEASURED CONSTANTS 

pf = fuel density, g/ml = 0.838 

Po = oil density, g/ml = 0.826 (15W-40) or 0.841 (SAE-30) 

(H/C)f = hydrogen to carbon ratio of fuel = 1.8 (light diesel) 

(H/C)o = hydrogen to carbon ratio of oil = 1.88 (15W-40) or 1.78 (SAE-30) 
LFE = Laminar flow element rating, CFM/in water 

E.lb CALCULATED CONSTANTS 

Kf = water formation constant of fuel 
= [g H 2 0/g H] x [g H/gfuel] 

= [(2 + 1 6)/2] x [( H/C )f/( 1 2 + (H/C)f)] 

= 1.174 

K 0 = water formation constant of oil 
= 1.219 (15W-40) or 1.163 (SAE-30) 

E.lc MEASURED VARIABLES 

Vf = volume flow rate of the fuel, ml/sec 

dP = differential pressure across laminar air flow element, inches water 

YY = amount of water condensed, ml 

SAcw' = specific activity of condensed water, dpm/ml 

A s = total activity of soot sample 

SAo' = specific activity of oil/kerosene mixture, dpm/(.5 ml of oil/kero solution) 
atm = barometric pressure, mm of mercury 



T a = air temperature, °C 
H = relative humidity, % 
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E.ld CALCULATED VARIABLES 

h s = specific humidity, g H20/g air 

= determined from psychrometric chart using H and T a 
Rf = wt rate of fuel, g/sec 
= Vf x pf 

(i c = viscosity coefficient 

= 1 - [.0025714 x (T a - 20)] (note — accurate for 15 <T a < 35) 
p a = density of air, g/m 3 

= (atm) x (133.32 Fa/mm Hg) x (28.962 g/mol) x (mol °K/8.3143 Nm) x 
[1/(273 + T a )°K] (note — ideal gas law) 

R a = wt rate of air, g/sec 

= (LFE) x (dP) x (.02832 m 3 /ft 3 ) x (min/60 sec) x [(atni)/29.92 in Hg] x 

(in/25.4 mm) x { 530 °R/[460 + (1 ,8T a + 32)]°R ) x pc x p a [8] 
SAcw = total specific activity of exhaust sample, dpm/g 
= SAcw' + (As/W) 

SAo = total oil consumption, g/hr 

= [RfKf + R a h s ]/[(SA C w/SA 0 ) - K 0 ] x (3600 sec/hr) 

E.2 OIL FLOW RATES (FROM FILM THICKNESS TRACES) 

RPM = revolutions per minute 

V 0 p = average oil volume difference, for a given region, between the 
power exchange strokes, mm 3 /cycle 
= calculated using an existing computer program which was modified 



for the Kubota (Appendix F) 
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V 0 j» = average oil volume difference, for a given region, between the 
gas exchange strokes, mm 3 /cycle 

= calculated using an existing computer program which was modified 
for the Kubota (Appendix F) 

V 0 = total avg. volume difference between the up and down strokes, mm3 
= ^op + 

Kop = rate of oil flow up the liner based on film thicknesses during the 
power exchange strokes, g/hr 
= (V 0 p) x (cnrVlOOO mm 3 ) x (RPM) x (60 min/hr) x 
(cycle/2 revolutions) x (p u ) 

Ros = rate of oil flow up the liner based on film thicknesses during the 
gas exchange strokes, g/hr 

= (Vq^) x (cm 3 /1000 mm 3 ) x (RPM) x (60 min/hr) x 
(cycle/2 revolutions) x (p 0 ) 

R</= overall rate of oil flow up the liner based on Film thicknesses, g/hr 
= Rop + Rog 

All of the preceding formulas are incorporated into LOTUS 123 v2.0 spreadsheets. 
Examples of these spreadsheets are included as the following six pages. The first two 
pages are for oil consumption calculations and the last four are for oil rates based on the 



film thickness data. 
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E.3 LOTUS 123 v2.01 SPREADSHEETS 



OIL CONSUMPTION CALCULATIONS: 

UNIFORM 



INPUT CALCULATIONS 



oil H/C ratio 1.78 

( SAE-30 = 1.78) 

( 15W-40 = 1.88) 

fuel H/C ratio 1.8 

activ of oil smpl 57052 

(dpm/.5 ml soln) 

oil density 0.841 

( SAE-30= .841) 

( 15W-40= . 826 ) 

Barometric Pressure 777 

(mm mercury) 

humidity coeff 0.01 

( lb wtr/lb air ) 

air temperature 28 

(Celsius) 



/ Oil formation constant - 

/ 

/ 

/ 

/ Fuel formation constant 
/ 

/ 

/ activ cone oil 

/ ( dpm/g ) 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 



-1 . 162554 



-1 . 173913 
-6919505 . 
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1 



UNIQUE 



INPUT / CALCULATIONS 

/ 

Load (lbs) 7.6 / TOTAL oil consumption 

/ (g/hr) 

Speed (rpm) 15 00 / 

/ 

/ 

fuel rate - C.2C33 / fuel rate 

(ml/sec) / (g/sec) 



air pressure diff 

(inches water) 


1.275 


activ of soot smpl 

{dpm/1 ml soln) 


JO? 


activ cor.c wtr sm.pl 

(dpm/ml) or (cpm/g) 


-i —j 


amount wtr collected 
(mi; 


2 . o 



/ 

/ air rate 
/ (g/sec) 

/ 



/ 


equivalence 


ratio 


/ 

/ 


(F/A)s = 0. 


0693 


/ 






/ 


total activ 


of soot smpl - 


/ 

/ 

/ 


( dpm ) 




/ 


total cor.c. 


of sample 


/ 


(dpm/ml) or 


(dpra/g) 



/ 

/ 

/ 

/ 

/ oil cor. sump (burned) 
/ oil con sump (ur.burn) 
/ Total 



1 . 42501S 

0 . 1 7 07 34 
3 . E 29377 
0 . 698160 
63 2 

- ■ 3 9 3 , o 2 

1.393132 
3 . 026734 
1.424916 
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OIL FLOWS UP THE LINER: 





PISTON 






AVG VOL 


OIL FLOW 


AVG VOL OIL FLOW 




REGIONS 






PWR STRK 


PWR STRK 


GS XCHNG GS XCHNG 


SAE-30 








( cu mm) 


( g/hr) 


(cu mm) 


(g/hr) 


1500 . 00 


I 


1/2 crwn lnd 


(X2) 


0.04 


1.51 


0 . 00 


0 . 00 




II 


top ring 




0 . 01 


0.38 


-0 .08 


-3.03 




III 


1 ring width 


below 1.11 


42.01 


-0 . 23 


-8.70 




IV 


2nd land 




2.30 


87 . 04 


-0.52 


-19 . 68 




V 


2nd ring 




-0.02 


-0.64 


-0 .11 


-4 .16 






PISTON 




AVG VOL 


OIL FLOW 


AVG VOL OIL FLOW 






REGIONS 




PWR STRK 


PWR STRK 


GS XCHNG GS XCHNG 


SAE-30 








( cu mm) 


(g/hr) 


( cu mm) 


(g/hr) 


3000 . 00 


I 


1/2 crwn lnd 


(x2) 


-0.34 


-12 . 87 


-0 . 12 


-4 .54 




II 


top ring 




0 . 08 


3 . 03 


-0 .21 


-7 .95 




III 


1 ring width 


below 0.74 


28.01 


-0 . 23 


-8.70 




IV 


2nd land 




1 .43 


54 . 12 


-0.45 


-17 . 03 




V 


2nd ring 




0.32 


12 . 11 


-0 .19 


-7 .19 







PISTON 


AVG VOL 


OIL FLOW 


AVG VOL 


OIL FLOW 






REGIONS 


PWR STRK 


PWR STRK 


GS XCHNG 


GS XCHNG 


15W-40 






(cu mm) 


( g/hr) 


(cu mm) 


(g/hr) 


1500 . 00 


I 


1/2 crwn lnd 


(x2) -0.04 


-1.51 


0 . 08 


3.03 




II 


top ring 


0.03 


1 . 14 


-0 . 04 


-1.51 




III 


1 ring width 


below 0.63 


23.84 


0.05 


1.89 




IV 


2nd land 


1 . 17 


44 . 28 


-0 . 08 


-3 . 03 




V 


2nd ring 


0.16 


6.06 


-0.01 


-0.38 






PISTON 


AVG VOL 


OIL FLOW 


AVG VOL 


OIL FLOW 






REGIONS 


PWR STRK 


PWR STRK 


GS XCHNG 


GS XCHNG 


15W-40 






( cu mm) 


( g/hr) 


(cu mm) 


(g/hr) 


3000 . 00 


I 


1/2 crwn lnd 


( x2 ) -0.22 


-8.33 


-0 . 16 


-6 . 06 




II 


top ring 


0 .10 


3.78 


-0 . 20 


-7 . 57 




III 


1 ring width 


below 0.46 


17 .41 


-0.38 


-14.38 




IV 


2nd land 


0 . 91 


34 . 44 


-0 .66 


-24 . 98 




V 


2nd ring 


0.09 


3 .41 


-0 .35 


-13 .25 
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OIL FLOWS UP THE LINER: 





PISTON 






AVG VOL 


OIL FLOW 


SAE-30 


REGIONS 






ALL STRK 


ALL STRK 








(cu mm) 


(g/hr) 


1500 . 00 


I 


1/2 crwn lnd 


(x2) 


0.04 


1.51 




II 


top ring 




-0 . 07 


-2.65 




III 


1 ring width 


below 0.88 


33.30 




IV 


2nd land 




1 . 78 


67 . 36 




V 


2nd ring 




-0.13 


-4.81 






PISTON 




AVG VOL 


OIL FLOW 






REGIONS 




ALL STRK 


ALL STRK 


SAE-30 








(cu mm) 


(g/hr) 


3000.00 


I 


1/2 crwn lnd 


(x2) 


-0.46 


- 17.41 




II 


top ring 




-0 .13 


-4 .92 




III 


1 ring width 


below 0.51 


19.30 




IV 


2nd land 




0.98 


37.09 




V 


2nd ring 




0 . 13 


4.92 






PISTON 




AVG VOL 


OIL FLOW 






REGIONS 




ALL STRK 


ALL STRK 


15W-40 








( cu mm) 


(g/hr) 


1500 . 00 


I 


1/2 crwn lnd 


(x2) 


0 . 04 


1 .51 




II 


top ring 




-0 . 01 


-0 . 38 




III 


1 ring width 


below 0.68 


25.73 




IV 


2nd land 




1.09 


41 . 25 




V 


2nd ring 




0.15 


5 . 68 






PISTON 




AVG VOL 


OIL FLOW 






REGIONS 




ALL STRK 


ALL STRK 


15W-40 








(cu mm) 


( g/hr) 


3000 . 00 


I 


1/2 crwn lnd 


(x2) 


-0.38 


-14 . 38 




II 


top ring 




-0 .10 


-3 .78 




III 


1 ring width 


below 0.08 


3 . 03 




IV 


2nd land 




0 . 25 


9.46 




V 


2nd ring 




-0 . 26 


-9.84 
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APPENDIX F DETERMINATION OF VOLUME DIFFERENCES 
F.l DETERMINATION OF REGIONS 

The program utilized uses a Crank-Nicolson technique to integrate the difference in 
areas under a film trace for the upstroke and downstroke pairs of specific regions. 
Subsequently, these area differences are swept around the circumference of the bore to 
determine the volume differences. 

Each inputted region is calculated based on using the "zero-point" crank angle as a 
reference position. This position, initially determined to equal -81 .6°, is defined as the 
angle at which the top of the piston passes the quartz window on an upstroke. The 
following relation is then used to determine the distance, s, in cm, from the crank axis to 
the piston pin axis at this point: 

s = a x cos(0) + (l 2 - (a 2 x sin 2 !©)) 1 / 2 

where: 1 = connecting rod length = 1 1 .01 cm 
a = crank radius = 3.5 cm 
9 = crank angle 

Next, a new "zero-point" crank angle is obtained by fitting a model of the piston to a 
printout of the film trace. Since piston geometry is known, this step allows for an estimate 
in the error of the original "zero-point". This error distance is subtracted (or added) to the 
original s and a new "zero-point' 'crank angle is determined by iterating. Specific regions 
can then be determined in a similar manner. 

For example: 

a) the initial "zero-point" is -81.6°. Therefore, s 0 = 10.96 cm 

b) it is determined from the trace that the s is 0.05 cm off. Therefore, 
the new s 0 = 1 0.96 - 0.05 = 1 0.91 cm 

c) by iterating, the new "zero-point" is determined to equal -82.5° 

d) the crown land is 0.95 cm long 

e) therefore, the new s = S] = s 0 + d = 10.94 + 0.95 = 1 1.89 cm 
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d) now iterate: 

0 = -67.7° implies St = 1 1 .85 
0 = -67.9° implies Sj = 11.84 
0 = -67.5° implies Sj = 11.86 (correct one) 

e) the region relating to the crown land would then be -82.5° to -67.5° 

f) the top ring region is determined in a similar manner except now 
s 0 = 11.91 cm 

Regions n, IV, and V consist of the top ring, second land, and second ring, 
respectively. They are used to reflect the film behavior of the upper portions of the piston. 
The crown land is not included based on the reasoning presented in section 6.2. Region I 
is determined based on the film on the liner as three centimeters of the crown land pass the 
window. The assumption is that the thickness on the liner at this point is the same up to 
TRR (top ring reversal). Finally, region III, one ring width below the top ring, is used in 
an attempt to substantiate the claim that oil is drawn up and around the top ring during 
operation. 
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C 

C 

c * • * 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

C • • • 

c 

c 

c 

c 

c 

c 

c 

c ••• 

c 

c 

c 

c 

c 

c 



PROGRAM MASSBAL 

Program to compute up / down stroke oil volume differences in 
distinct specified regions along the piston for multiple Sloan 
data files and write them to an output file that can be 
converted to a bar graph format. 

This program requires as input any number of standard Sloan data 
files containing an array of position vs film thickness. The program 
will average the volume of oil in any number of user specified 
regions along the pistan. 

The program has been adapted to analyze either McElwee 
[where IRDATE(3)-89 ] ar Bliven flar.* data sets. Commands 
specific to Lux [where IRDATE(3)— 88 ] data sets have been 
neutralized by comment lines which are specified by "cc". 

[M. Bliven, 14 March 1990] 



Program structure: 

Main Level 1st Sublevel 2nd Sublevel 3rd Sublevel 



MASSBAL 




FSTDCTPT 

INTECRATEDF 

AVERAGE 






I/O 


Un i t 


3: 


INFILE Slaan 


data file 


Un i t 


5: 


Terminal 


I 




Un i t 


6: 


Terminal 


0 




UNIT 


10: 


CXJTFILE, 


BDC 


LOCATIONS 



i mpo r t an t variables 



AVALS 

CALOC 

NPMAX 

CALFCTR 

CAOFTST 

FI LM(2 . NPMAX) 

GEMNVAL 

GESD 

INFILEV 

INTRPM 

LWRB (10) 

MNVAL 

NBDCLOC 

NSTRCKE 

NUMFIL 

NUMREG 

OFFST 

PYrVNVAL 

PWSD 



Array af average volume for like revolutions 
Crank angle location 

parameter, max number af points to plat 
calibration factor far individual data set 
correction factor for shaft encoder offset 
Array of position vs film thickness 
Mean value of oil volume difference for gas 
exchange strokes 

Stondard deviation for oil volume difference 
for gas exchange strokes 
Input file vector 

Value for engine speed calculated from time 
difference between BDC pulses 
Lower bound for region i 

Mean value of oil volume for a region for like 
revolutions from the program AVERAGE 
Botton dead center location in crank angle 
Number of strokes recorded in data file 
Number of Sloan files to process 
Number of regions ta divide data into 
Instrument offset factor to initialize zero 
film thickness 

Mean value of oil volume difference between 
powered strokes 

Standard deviation or oil volume difference 
between powered strokes 



o o o o o o o 
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SIGVAL 



REGVAL 



Arroy contoining piston regions with the 
overoge volue of oil volume difference for 
I i ke revo I u t i ons 

Stondord deviotion of oil volume for like 
revolution from the progrom AVERAGE 



PARAMETER NPMAX - 190000 
PARAMETER NPCMAX - 60000 
PARAMETER NCMAX - 250 
PARAMETER NREGMAX - 15 



I mox / of points/ch 
I mox / of points/cycle 
I mox/ of cycles 



C 



INCLUDE ’ SLNCCM. FOR * 



I contains IMPLICIT INTEGER-4 



C 



C 



INTEGER-2 IHED1 (256) , IHED2(256) 
EQUIVALENCE (IHED1 .NRUN) , ( IHED2 , NS LOT) 

INTEGER-2 IDATA(NPMAX) 



C 



INTEGER NBDCLOC(NCMAX) .NCHOICE. NSTROKE( NCMAX) , 

$ NUMREV.PCH, ICH.Q. ICH1 ,CALOC(365) .NUMREG 

REAL F I LM ( 2 , NPCMAX ) , AVALS (2 , NCMAX) , VALUE(2). 

$ XOLD.XNErf, INTRPM(NCUAX) , CALFCTR , ROLOFF .CAPPT , 

$ XSCL(4) , INTAR (10000) , CUTAR ( 1 0000) .DUMC7, 

$ SUM , VAL, OFFST , WVAL, UHPCS . SIGVAL, S ICPOS , DUMPLOT , 

$ LSTROKE, CKRAD, CRDL , BORE. CAOFFST, LWRB(NREGUAX) , 

$ UPRB(NREGUAX) , PYr%*fVAL(NRECMAX ) , PWSD (NREGMAX ) , 

$ CEMNVAL(NREGUAX) , GESD(NREGMAX) , REG VAL (NREGMAX. NCMAX) 



CHARACTER INF I LEV(20) -40 , OUTFI LE- 40 , A- 1 , B- 1 , C- 1 , 
$ CALI NT -3 . STROKE- 4 .skip-3 
C 

C • • • initialize: 

C 

CALINT - ‘ON* 

ICH - 3 



WRITE (5,1000) CALINT 
READ (5,-, ERR - 10, END - 999) NCHOICE 
C 

C ••• NCHOICE - 1: Set dota files to read 
C 

IF (NCHOICE. EQ . 1 ) THEN 

WRITE (6 , * ( * ' How mony flor.- files? (up to 20): ’*,$)') 

REA0(5.-) NUMFIL 
DO 1-1, NUMFIL 

WRITE ( 6 , * ( ’ * Enter f i I enome / M ,I2, M : # *$) # )I 
READ(5 , 1 (A40) ’ ) INFILEV(I) 

ENDDO 

ENDIF 

C 

C ••• NCHOICE - 2: Set piston regions by cronk ongle 
C 

IF (NCHOICE. EQ.2) THEN 

WRITE (6,'(*’ How mony piston regions? (up to 16): '*,$)*) 
REA0(5,e) NUMREG 
DO 1-1, NUMREG 

WRITE (6,’(*’ Enter lower, then upper bound (CA deg) for 
$ region / M ,I2, M : *•$)•)! 



C 



C 

C - - - set up menu : 



C 

10 CONTINUE 



ISTAT - LIB$ERASHLPAGE(1 , 1) 




ooooooo ooo 
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REA0(5 , • ) IWRB(I), UPRB(I) 

ENDOO 
END IF 

••• NCHOICE — 3: Linearly interpolate between crank angles 
IF (NCHOICE. EQ. 3) THEN 

WRITE(6 . ’ ( IX’ ' Li near ly interpolate between crank angles? 




IF (C.EO. * Y * . OR . C . EG . ' y * ) CALINT-'ON' 
IF (C.EQ. # N * . OR . C . EQ . ’ n ’ ) CALINT-'OFF* 
ENDIF 



NCHOICE - 4: Chock entered information ••• 

Note: The running time of the program could be lengthy 
depending on the number of data files and regions specified. It is 
wise to check the Input data for accuracy prior to execution. 

IF (NCHOICE. EO. 4) THEN 

ISTAT - LIBS ERAS E_PAG E ( 1 , 1 ) I erase page 
DO I - 1 , NUMFI L 

WRITE ( 6 . ' ( * * Data File / 'M3," is: ‘\A40)’)I, 

$ INFILEV(I) 

ENDOO 

DO I - 1, NWREC 

WRITE (S,’(’* Crank angles bounding region f * , 12, 

$ ,B are: * ’ # F7. 2. * * to ’ * . F7 . 2) * ) I , LWRB(I), UPR0(I) 

ENODO 

WRITE (6.*(*' RTN to continue '*,$)*) 

READ(5 . ’ (A 1 ) * ) B 
ENOIF 
C 

C ••• NCHOICE - 5 Execute main program and write output. 

C 

IF (NCHOICE. EO. 5) THEN 
DO I - 1 . NUUFIL 

ISTAT - IIB$ERASE_P'4GE( 1 1 ) lerase page 
WRITE (5.1001) CAIINT, INFILEV(I) 

ILUN - 3 

OPEN(UNIT-ILUN,NAME-INFILEV( I ) .TYPE-’ OLD’ , ACCESS- *0 IRECT ’ 

$ .FOR)*-’ UN FORMATTED’) 

C 

C ••• and read header: 

C 

READ ( I LUN * 1 ) (IHEDI(L).L-I .256) 

READ ( I LUN ' 2) (IHED2(L).L-1.256) 

C 

I DU V - 1 

MPPC - IDUW*NBPC*256 
C 

C ••• load up the private common 
C 

C These values are user specified prior to data collection. 

C 

CAPPT - RUSER(26) 1 crank angles per point 

cm' 
cm' 
cm 
C 

C ••• Convert to mm 



LSTROKE - RUSER(27) 
BORE - RUSER(28) 

CRDL - RUSER(29) 



I engine stroke 
1 engine bore 
I connecting rod length 
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c 

BORE - BORE* 10. 

LSTROKE - LSTROKE* 1 0 . 

CRDL - CRDL* 1 0 
CKRAD - LSTROKE /2. 

C 

C Sot shaft oncodor error in variable CAOFFST 
C 

C The shaft encoder marks the crank angle at BOO. Depending on 

C certain errors the recorded BDC position is not necessarily at -180. 

C 

cl IF (IRDATE(3).EQ.88) THEN 

cc CAOFFST - -178.5 

cc ELSE 

if ( i rda t e (3) . eq . 89 . and . n run . gt . 99 . and . n run . I t . 1 06) then 
caof f st - -181 
else 

if ( i rda t e (3) . eq . 90) then 
caof f s t - -180.926 
e I se 

CAOFFST - -180 
end i f 
end i f 

cc END IF 

C 

C ••• Look for old BDC file. The BDC file contoins the specific 
C calibration foctor and the number of revolutions in the data set. 

C 

OLITFILE - * * 

ENCOOE(40, *(• ’ BQCLOC . * M3) \OUTFILE)NRUN 
C 

C ••• Open existing BDC file 
C 

CPEN(UNIT-10.NAUE-CUTFILE.TYPE-‘OLD* .ACCESS-* SEQUENT IAL * , 

$ FORU-’FORUATTED' ) 

READ( 10.*) CALFCTR. NUUREV 
DO J-1 .NOJAX 

READ( 10. • , ENO-15.ERR-5) NBDCLOC(J) . INTRPW(J) .NSTRCKE(J) 

ENDDO 

15 CONTINUE 
C 

C ••• Set region along piston to analyze 
C 

DO N - 1, 2 
DO J-1 ,NL%«EG 
DO L-1 , NUM3EV 

REGVAL(J.L) - 0.0 
ENDDO 
ENDDO 
M - 1 

DO K-N, NUWREV-1 . 2 
IFPT— N6DCLOC (K) 

C 

C ••• Set I FPT and NRPTS Set PCH for labels 
C 

NRPTS - NBDCL0C(K+1) - NBDCLOC(K) +1 I / of points to read 
IF ( NRPTS. GT.NPMAX) NRPTS - NPUAX 
C 

C •*• Linearly interpolate between crank angles 
C 

IF (CALINT. EQ. ‘OFF') GOTO 21 
I CHI - 1 
C 
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c 



c2 

cc 

C 



c3 

cc 



20 



21 

C 

C • • • 

c 

c 

c 1 
cc 
cc 
cc 
cc 
cc 

c 



c 



CALL FSTDGTPT ( I LUN , ICH1 , IFPT .NRPTS , I DATA) 

NUMCA - 0 

XOLD - CSCL(ICH1)*IDATA(1)+CBIA(ICH1) 

IF (NRUN. LT.47. and . i rdate(3).eq.88) XOLD-XOLD-5 . 0 
IF (NRUN. GT. 47. and. i rdot e(3) . eq . 88) X0LD-X0LD+5 . 0 
skip « ' no ’ 

DO L-2, NRPTS 

XNEW - CSCL( ICH1 ) • I DATA( L)+C8 I A( 1CH1 ) 

IF (NRUN. LT. 47. ond . I rdot 9 (3) . eq.88) XNEW-XNEW-5 . 0 
IF (NRUN. GT. 47. and. i rdot«(3) .«q.88) XNEW-XNEW+5 . 0 
IF (XNEW. GT. 2. 5. AND. XOLD. LT.2.5) THEN 
if (capp t . «q. .5. and . sk i p. «q . ’y as ' ) thsn 
skip - * no ’ 
go to 20 
and i f 

if (capp t . eq . .5 . and . ak i p . eq . ’ no ') skip-'yes* 
NUMCA-NUMCA-fl 
CA LOC(NUWCA)-L + IFPT 
ENDI F 
XOLC^XNEW 

ENODO I END OF L LCCP FOR CRANK ANGLE LOCATION 

CONTINUE 



Average first 100 data points of pmt signal ta 
estaolish zera affsat. .005 sets this ta 2 raicrans 



NOTE: This section adapted to analyze data collected 

by BLIVEN whose PVT signal due to instrument offset ond 
laser reflectances alane (ie. dry liner) was measured to be 
approximately 0.008 V. Averaging over first 100 data paints 
is by-passed and offset value is set ta 0.008. 

CALL FSTDGTPT (I LUN, ICH , IFPT .NRPTS , I DATA) 

if ( i rda t e (3) . eq . 89 ) then 

SIM -0. 

DO L ^300 , 400 

VAL - CSCL(ICH)*IDATA(L) + CBIA(ICH) 

SIM - SUM + VAL 
DLMCT - L— 299 

ENDDO I END OF L LOOP FCP OFFSET 

OFFST - SLM/DUUCT— .016 I zero pt offset for WcElwee 
else 

OFFST - .0 I zero pt offset for Bliven 

end i f 



Convert the data ond write them into plot-orray 

Note: CLKRTE is clock-rate in [kHz] and con be used to 
convert x— axis data time 



DUMP LOT- CAOFFST 
Ol 

DO L - 1 .NRPTS 

••• Convert time to crank angle by speed or position Interpolation 
IF (CALINT . EQ. ’OFF ’ ) THEN 

FILM(I.L) - CAOFFST + 260 . ‘FLOAT ( L— 1 )/FL0AT (NRPTS— 1 ) 

ELSE 
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c 

C *** If Cronk angle pulse is coincident with sample number, then 
C increment x value by one crank angle 

C 

IF ( (L+IFPT) . EQ.CALOC(Q) ) THEN 
FI LM( 1 , L) - CAOFFST + FLOAT (Q) 

DUMPLOT-FI LM( 1 , L) 

00+1 

ELSE 

C 

C *** Otherwise divide fractions of cronk angle evenly between pulses 
C 

IF (0.E0.1) THEN 

FI LM( 1 , L) - CAOFFST+FLOAT ( L— 1 )/FLOAT (CALOC(Q)-IFPT) 

ELSE 

FI LM( 1 . L)-0UMPL0T+1 . 0/FLOAT( (CALOC(Q)-CALOC (0-1 ) ) ) 
DUMPLOT-F I LU( 1 , L) 

END IF 
END IF 
ENDIF 

FI LW(2 , L)-(CSCL( ICH) * IDATA(L)+C9IA( ICH)-OFFST)/CALFCTR 

ccc 

ccc wri te(6, * ( * * FI LM(2 , L) - ‘\M0.5)*) f I lm(2. I ) 

ccc wr i te(6, ’ ( 2 i 3 . i6,2f8.3) ')n t k. I . f i I m ( 1 . I) , f i Im(2. I ) 

ccc 

ENDDO I END OF L LOOP POP FILM ARRAY 
C 

C ••• Coll integration routine for eoch region. Program INTEGRATEDF 
C will calculate the volume of oil In o specified region for one 

C revolution. 

C 

DO J-l . NUMREG 

CALL INTEGRATEDF ( FI LM . LWR9( J ) .UPRB (J ) .VALUE .NRPTS . CROL. 

S CKRAO) 

REGVAL(J.M) - VALUE(2)*BORE*3. 14159 

ccc 

ccc write(6,’(4i3,3e10.4)')m,n ( k,j,volue(1),value(2),regvol(j,ffl) 
ccc 

ENDDO I END J LOOP FOR REGIONS 
AVALS( 1 .M) - VALUE(I) 

M-a+i 

ENDDO I END OF K LOOP FOR ALL REVOLUTIONS OF 1 TYPE 
C 

C ••• Coll overoge routine for like revolutions of one region. Program 
C AVERAGE will colculote the overoge volume for the like revolutions. 

C 

DO J-1 .NUMREG 

DO L - 1 . N Uk«EV 
AVALS (2. L) - 0. 

ENDDO 
DO L-1 .U 

AVALS(2.L) - REGVAL(J.L) 

ENDDO 

CALL AVERAGE( AVALS . MNVAL . MNPOS . S IGVAL . S IGPOS ) 

IF (NSTROKE(N) . EQ. 1 ) THEN 
PWMNVAL(J) - MNVAL 
PYfSD(J) - SICVAL 
ELSE 

CEMNVAL(J) - VWVAL 
GESD(J) - SIGVAL 
ENDIF 

ENDDO I END OF J LOOP (DIFFERENT REGIONS) 

ENDDO I END OF N LOOP (DIFFERENT REVOLUTION TYPES) 
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C 

c *** Write output to o file with format mossbol./ file 
C 

ENCODE (40 , ' ( ” [ 0 LIVEN. DATA. MASS ]MASS8Al. ” , 13) ' , OUT FILE) NRUN 
OPEN (UN IT- 10 . FILE-OUTFILE,TYPE-'NEW’ .ACCESS- 'SEQUENT IAL ' . 

$ FORM- 1 FORMATTED*) 
c 

write (10. '(lx” ”)') 

write (I0,’(28x”| — — — — I ”)’) 

write (10, * (28x * * j MASS8AL RESULTS | ")') 
write (10,’(28x”l ■ | ")') 

write ( 1 0 , ' ( 1 x ' * ' ' ) ’ ) 

write ( 1 0 . ’ ( 2Sx ” Da t o Set: ’ ’ , al5) ' ) inf i I ev( \ ) 
write (10. •(lx” ' ')') 
write (10. ’(lx** ")') 
c 

WRITE(10, ' (5X * * POWER STROKES GAS 

$ EXCHANGE ” ) ' ) 

wri te(l0, ' ( 5 x ' * 

$ '•)’) 

WRITE( 10,' (5X ’ ’REG FRM TO MEAN STD DEV MEAN 

$ STD DEV”)') 

WRITE( 1 0 , ’ (5X * * f (deg) (deg) (cu mm) (cu mm) (cu mm) 

$ (cu mm) ” ) • ) 

wr i t e( 1 0 , ' (5x ” 

$ ••)•) 

DO L -1 . NtAtfEG 

write (10, ’(lx” ” ) ’ ) 

WRITE (10. ’(5X. 12, 1X,F7.2.2X,F7.2.4(3X.F7.2))')L. LWR0(L), 

$ UPR0( L) . PWMWAL(L). PWSD(L). GE)^VAL(L), GESD(L) 

ENCDO 

write (10, ’(lx* * ”)') 
write (10, ’(lx” ” ) * ) 

write ( 10, ’ ( 1 x” ” ) ' ) 
write ( 1 0 , ' ( 1 x * ' ” ) ’ ) 

write ( 1 0 , * ( 1 x ” NOTE: If negative crank ong I e regions specified. 

$ then POSITIVE mean values indicate ”)*) 

write (10, ’(lx” that UPstroke volumes are GREATER than 

$ DCWNstroke volumes.'')') 

write (10. ’(lx” Standard deviotions are always positive.”)') 

END DO I END OF I LOOP (MULTIPLE FILE) 



c 


END IF 1 END OF NCHCICE 


5 IF LOOP 


c • • • 

c 


Final Branch 




5 


CONTINUE 




c 


IF (NCHOICE.NE. 0) GOTO 


10 


999 

C 

c 


CONTINUE 




c * • * 

c 


Formats: 




1000 


FORMAT (/' ••• OIL FILM 


STATISTICAL ANALYSIS •**',/// 



ST 10, ’ ( 1 ) Enter files to analyze. ' ,/, 

$T10,'(2) Set piston regions. * ,/, 

$T10,'(3) Llneor Cronkongle Int e rpo I at i on : ’ , 1 x , A3 ,/, 
$T10,'(4) Check Information for 1 k 2. ' ,/, 

$710, '(5) Compute averages ond write to massbal file.',/, 
$T 1 0 , ' (0) or CTRLtZ for exit’///. 

$ ' Enter cho ice: ’ ,$) 

1001 FORMAT (/’ ••• OIL FILM STATISTICAL ANALYSIS •••*.///, 
$T10.'(1) Enter files to onolyze. ’ ,/, 
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$T10,'(2) Set piston regions. " ,/. 

$T1 0 , * (3) Linear Crankangle Interpo I at i on : 1 , 1 x ,A3 ,/, 
$T 1 0 , ' (4) Check information for 1 k 2. ' ,/» 

$T10,'(5) Processing file: *^30, * Be patientl',/, 

$T 1 0 , ’ (0) or CTRLtZ for exit*///. 

$’ Enter choice: *,$) 

STOP 

C 

END 
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F . 3 "INTEGRATEDF.FOR" LISTING 
C 

SUBROUTINE INTEGRATEDf ( FI LM , LED , UBD , VALUE , NRPTS . CRDL , CKRAD ) 

C 

C ••• This suboutine integrotes the volue in Film(2,I) using the Cronk- 
C Nicholsen method. It reods the step size from the first row of 

C Film. The ronge is user-specified with lower ond upper bounds. 

C The difference between oreos of successive strokes is found and 

C returned to the driver program for overaging. 

C 

C ••• DIMENSION AND DECLARE VARIABLES 
C 

PARAMETER NPMAX - 250000 ! max / of points/ch 

PARAMETER NPCUAX - 60000 I mox f of points/cycle 

PARAMETER NCMAX - 250 1 mo xf of cycles 

C 

REAL FI LM(2 , NPCMAX) , LSD , UBD , VALUE(2) , AREA, POSA, POS3. 

$ CRDL , CKRAD , DMBD 
INTEGER I , J , K , NRPTS 
C 

C ••• SET INITIAL CONDITIONS BEFORE INTEGRATING 
C 

J - 1 
6 K - 1 

DO 1-1. NRPTS 

IF (LBD.CE. FI LM( 1 .1)) THEN 
K-I 

ccc 

ccc write(6,’(2M0.3)*)lbd,film(l,k) 

ccc 

ELSE 
GOTO 3 
END IF 
ENCDO 

3 CONTINUE 

C 

C ••• Perform integration 
C 

AREA - 0. 

10 IF (FI LW ( 1 , K+1 ) . LE. UBD) THEN 

POSA - CKRA0*C0S(FILM(1 ,K)/360. *20. 141593)+ 

$SCRT (CRDL* *2 - (CKRAD* *2) • (SIN(FI LM( 1 ,K)/360. *2*3. 1 41593) ) **2) 

POSB - CKRAD*C0S(FILM(1 , K+1 )/360. *2*3. 14159)+ 

$SCRT (CRDL* • 2 - (CKRAD* *2) • (S IN(FI LM( 1 ,K+1 )/360 . *2*3 . 1 41 59) ) • *2) 

C 

C ••• Divide by 1000. to convert microns to mm 
C 

AREA-AREA+( (FILM(2,K)+FILU(2,K+l))/1 000 . ) • 

$ A8S ( (PQS3— P0SA)/2 . ) 

K - K+1 

ccc 

ccc wri te(6, ' (3e10.3) *)po3a,posb,areo 

ccc 

GOTO 10 
ENDIF 

ccc 

CCC write(6, *(2f10.3,e10.3) , )fi lm( 1 , k) , f i lm(2 . k) . area 

ccc 

C 

C ••* Mirror bounds ond find the oreo of the adjacent downstroke, 

C then subtroct the values ond return. 

C 

IF (J.EQ.1) THEN 
DM80 - LSD 
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LBO - -UBD 
USD - — 0M8D 
VALUE(I) - AREA 
J - 2 
GOTO 6 
END IF 

VALUE(I) - VALUE(I) - AREA 
VALUE(2) - VALUE(I) 

•** Return bounds to original values 

DMBD - LBO 
LBD - —UBD 
UBD - -OWED 

• That * s the end 

RETURN 

END 



o o o o o o 
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F . 4 "AVERAGE. FOR” LISTING 

C 

SUBROUTINE AVERAG E( AVALS .MNVAL .MNPOS t S IGVAL , SIGPOS) 

*** SUBROUTINE COMPUTES AVERAGE VALUES AND STANDARD DEVIATION 
FOR THE FEATURE OF INTEREST N PASSED IN MATRIX * AVALS * 

*** DIMENSION AND DECLARE VARIABLES 

PARAMETER NPMAX - 250000 1 max / of paints/ch 

PARAMETER NPCMAX - 60000 ! max / of points/cycle 

PARAMETER NCMAX - 250 l max/ of cycles 

C 

REAL AVALS(2. NCMAX) .MNVAL . MNPOS , S IGVAL . S IGPOS 
$ XSUU, SUMPOS , XVAL , CMCT , XDEV. POSDEV, SDSUM, SDPSUM 
INTEGER I.J.K 
1-1 

XSUU-0. 

SUMPOS— 0 . 

10 IF (AVALS (2 , I ) .NE. 0. 0) THEN 

XVAL - AVALS (2 , I ) 

XSUM - XSUU + XVAL 
SUMPOS - SUMPOS + AVALS (1.1) 

DMCT - I 
I-I + 1 
GOTO 10 
END IF 

MNVAL - XSLM/CMCT 
MNPOS - SL^CS/CMCT 
C 

1-1 

SDSLM-0. 

SCPSUM-0. 

12 IF(AVALS(2, I ) .NE.0.0) THEN 

XDEV - (AVALS (2.1) -**WAL)-«2 
POSDEV - (AVALS (1,1) - WMPOS)««2 
SDSLM - SDSUM + XDEV 
SDPSUM - SCPSUW + POSDEV 
CMCT - I 
I-I + 1 
GOTO 12 
END IF 

SIGVAL-SCRT(SDSLM/CMCT) 

S IGPOS— SORT (SDPSIM/DUCT ) 

RETURN 



C 



END 
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Figure 1 . Engine Instrumentation/Measurement Equipment 
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Figure 2. Piston and Rings of the Kubota EA300N 
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Measured Measured 




Figure 3. Schematic of Sampling System 
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Figure 4. Procedural Flow Chart 



oil consumption (g/h r) 
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OIL CONSUMPTION USING SAE-30 
(full load @ 1500 and 3000 rpm) 

3.0 -| 



2.6 - 



2 . 2 - 



1.8 - 



1.4 - 



1.0 - 



0.6 

1000 1500 2000 2500 3000 

RPM 



□ 



0 

□ 



D 



1 » I 



3500 



Figure 5 



oil consumption (g/hr) 
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OIL CONSUMPTION USING 15W-40 
(full load <§> 1500 and 3000 rpm) 




Figure 6 



oil consumption (g/hr) 
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GRAPH OVERLAY OF SAE-30 AND 15W-40 
OIL CONSUMPTION 




1 5W-40 
SAE-30 



Figure 7 



oil consumption (g/hr) 
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COMPARISON OF THE OIL CONSUMPTIONS 
USING SAE-30 AND 15W-40 



■ 15W-40 
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Figure 8 



oil consumption (g/hr) 
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VOLATILITY OF NEW SAE-30 
(full load @ 1500 rpm) 




Figure 9 



FILM THICKNESS (mcrons) 
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STROKE : COMP 
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Figure 10 



FILM THICKNESS (mcrons) 
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STROKE :EXP 
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Figure 11 



FILM THICKNESS (mcrons) 
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STROKE :EXH 
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Figure 12 



FILM THICKNESS (microns) 
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STROKE : I NT 
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Figure 13 
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Figure 14. Areas /Volumes of possible oil 
consumption mechanisms 
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LABELED REGIONS OF COMP STROKE: 
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FILM THICKNESS (mcrons) 
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LABELED REGIONS OF COMP STROKE: 
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Figure 16 



OIL FLOW UP THE LINER (g/hr) 
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OIL FLOW FOR THE POWER EXCHANGE STROKES 
USING SAE-30 AT 1500 AND 3000 RPM 
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Figure 17 
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OIL FLOW FOR THE POWER EXCHANGE STROKES 
USING 15W-40 AT 1500 AND 3000 RPM 
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OIL FLOW UP THE LINER (g/hr) 
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OIL FLOW FOR THE GAS EXCHANGE STROKES 
USING SAE-30 AT 1500 AND 3000 RPM 
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OIL FLOW UP THE LINER (g/hr) 
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OIL FLOW FROM THE GAS EXCHANGE STROKES 
USING 15W-40 AT 1500 AND 3000 RPM 




Figure 20 
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OIL FLOW FOR ALL STROKES 
USING SAE-30 AT 1500 AND 3000 RPM 
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Figure 21 
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OIL FLOW FOR ALL STROKES 
USING 15W-4G AT 1500 AND 3000 RPM 

100 



80 

k. 

x: 

3 

a: 60 

LU 

Z 



x 40 

a 

i 20 



O 0 



-20 



■ 1500 RPM 15W-40 

□ 3000 RPM 15W-40 




oil cons. I II III IV V 



PISTON REGION 



Figure 22 



Ro/Ro' (Ro’= respective value at 1500 rpm) 
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NON-DIMENSIONALIZED OIL CONSUMPTION 
AND REGION II VALUES FOR POWER 
EXCHANGE STROKES USING SAE-30 




Figure 23 
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NON-DIMENSIONALIZED OIL CONSUMPTION 
AND REGION II VALUES FOR POWER 
EXCHANGE STROKES USING 15W-40 
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OIL FLOW RATE UP THE LINER (g/hr) 
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OIL FLOW RATES OF THE UPPER 
HALF OF THE PISTON 
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TEMPERATURE (Celsius) 
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RPM vs. UPPER LINER TEMPERATURES 




Figure 26 
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VISCOSITIES OF SAE-30 AND 15W-40 

vs. RPM 
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Figure 28 
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VISCOSITIES vs. OIL CONSUMPTIONS 




Figure 29 
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